Impact of the Dzyaloshinskii-Moriya interaction in strongly correlated itinerant 

systems 
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Spin-only approaches to anisotropic effects in strongly interacting materials are often insufficient 
for systems close to the Mott regime. Within a model context, here the consequences of the low- 
symmetry relevant Dzyaloshinskii-Moriya (DM) interaction are studied for strongly correlated, but 
overall itinerant, systems. Namely, we investigate the Hubbard bilayer model supplemented by a DM 
term at half filling and in the hole-doped regime. As an add-on, further results for the two- impurity 
Anderson model with DM interaction are also provided. The model Hamiltonians are treated by 
means of the rotational invariant slave boson technique at saddle point within a (cellular) cluster 
approach. Already small values of the anisotropic interaction prove to have a strong influence on the 
phases and correlation functions with increasing U. An intriguing metallic spin-flop phase is found 
in the doped bilayer model and a reduction of the RKKY exchange in the two-impurity model. 

PACS numbers: 71.27.+a, 71.23.An, 75.30.Hx, 75.20.Hr 



I. INTRODUCTION 



The effect of anisotropic magnetic exchange on the 
atomistic level has been recently brought back to a cen- 
tre of interest in condensed matter physics due to its 
intriguing importance in e.g. the search for multifcr- 
roic materials^, the understanding of complex metal- 
lic magnetic surface structures 3 - or the phenomenology of 
topological insulators^ A hallmark step in this research 
topic has already been performed some fifty years ago 
by Dzyaloshinskii^ and Moriya^ who derived an effective 
spin-spin interaction term from the spin-orbit coupling in 
low-symmetry cases, the so-called Dzyaloshinskii-Moriya 
(DM) interaction. The DM term or more generic mag- 
netic anisotropics are nowadays believed to play further- 
more a prominent role in many strongly correlated ma- 
terials. However nearly exclusively, theoretical studies in 
this context were in the past bound to pure spin models 
without itinerancy, leaving the impact of charge fluctua- 
tions aside. Yet the latter are surely important, e.g. close 
to the Mott-critical regime of the metal-insulator transi- 
tion. Allying the Hubbard model with spin-orbit terms 
has just recently gained rising interest^ 

In the present work we aim at a minimal modeling of 
the influence of the DM interaction in the strongly corre- 
lated metallic regime. There are many specific materials 
problems motivating such a case study, namely the com- 
plex magnetic behavior of doped cuprate systems 
manganiteS ) 12 ' 13 and mono-oxides^ as well as anisotropic 
magnetic effects close to the metal-insulator transition 
in low-dimensional organic compounds^ or in the con- 
text of transition-metal impurities on metallic surfaces^ 
While standard direct and indirect exchange processes 
favor collinear alignment of the local spins generated in 
the strongly correlated metallic regime, the DM interac- 
tion tends to align the spins in a perpendicular fashion. 
Thus the competition between the former conventional 
exchange processes and the DM interaction within an 
itinerant system shall give rise to nontrivial physics re- 



sulting in sophisticated spin arrangements/orderings. 

To keep things simple and to build up on a somewhat 
canonical approach, we rely on two basic models, namely 
the bilayer model of two coupled single-band Hubbard 
planes^— and the two-impurity Anderson model^ - — 
The former Hamiltonian allows for a DM coupling be- 
tween two lattice planes in the thermodynamic limit, 
whereas the latter one provides the possiblity to study 
the DM term within a local perspective via interact- 
ing impurities coupled to the same bath. Both setups 
render it possible to investigate nearest-neighbor (NN) 
correlation functions between sites in an itinerant back- 
ground. Of course, such modelings are not sufficient 
to grasp the very details of the above named materials 
problems, yet it will be shown that the computed phe- 
nomenology is far from trivial and may apply to generic 
realistic phenomena. One key focus in the context of 
the Hubbard bilayer lattice is thereby on the competition 
between the antifcrromagnctic (AFM) tendencies driven 
by direct exchange and the DM term within the metal- 
lic state. It will become clear that already rather small 
values of the DM integral may have a significant influ- 
ence on the magnetic ordering tendencies in the larger 
Hubbard U range, i.e., the AFM state is rather sensi- 
tive to only minor DM perturbations. A rich phase dia- 
gram results from the interplay of kinetic energy, onsite 
Coulomb and DM interaction. The latter also has im- 
portant consequences in the two-impurity model, where 
its favor for perpendicular spin arrangement severely af- 
fects the local-limit competition between singlet-forming 
Kondo-screening and triplet-forming Ruderman-Kittcl- 
Kasuya-Yosida (RKKY) interaction. 

In the following we define the model Hamiltonians as 
well as our mean- field approach in section [XT] The results 
for the Hubbard bilayer at half filling and in the hole- 
doped case are discussed in section IIIII Some basic ob- 
servations retrieved from the studies on the two-impurity 
Anderson model with DM interaction will be presented 
in section II VI 



2 



II. HAMILTONIANS AND THEORETICAL 
APPROACH 

The first problem addressed here consists of two cou- 
pled two-dimensional infinite square-lattice planes with 
one orbital per site each facing an on-site Coulomb repul- 
sion U (see Fig.[T]). In both planes the electron dispersion 
is defined by identical simple NN hopping t. The inter- 
plane coupling is realized via a perpendicular hopping t± 
as well as a DM interaction mediated by the vector in- 
tegral D. The model Hamiltonian is accordingly written 
as 

, 0<T , ia 

+ ^D.(S i; xS 2i ) , (1) 

ai i 

where c^J a creates/annihilates an electron in layer a=l,2 
at lattice site i with spin projection cr=t, \» The u=x, y, z 
component of the spin operator at each site i of an indi- 
vidual layer a is provided by S£/=l/2 icr r„}, c aia , with 
the Pauli matrices t^. In general, the vector interaction 
D is defined perpendicular to the bond between the in- 
volved lattice sites££ Since otherwise there is a freedom 
of choice for the explicit direction, we pick D to point 
along the y axis, i.e. T)=De y . Note that the DM in- 
teraction may only occur if the inversion symmetry is 
broken. To facilitate this in the present case, one could 
e.g. think of an inter-layer coupling originally established 
via oxygen with an angle deviating from 180°. 

In the smaller second part of this paper, we take the 
opportunity to also briefly discuss the well-known two- 
impurity Anderson model (TIAM) supplemented by a 
DM interaction between the impurities. We write that 




FIG. 1. (Color online) Bilayer model with DM interaction. 
The vertical blue line is a representant of the inter-layer hop- 
ping t±, which of course is applied at every lattice point, and 
the ellipse marks the two-site cluster. The DM integral vector 
D is chosen to point in y direction. 



model in the form 

#TIAM ^ c la C ka + £ d^2 nia 

ktr i(J 

+ VY^tLdi* + h - C + *12 X^L^a + h.C.) 

+ l/^n ltf!li + D-(S 1 xS 2 ) , (2) 

i 

with the impurity-electron operators df^J (i=l,2), the 
impurity-level energy Ed and the impurity-impurity hop- 
ping t\2- The bath has associated operators cfj and a 
dispersion e^. The impurity-bath coupling is denoted by 
V and the Hubbard U is located on the impurities with 
m a =d\ a d ia . In the present approach the bath is treated 
explicitly through a three-dimensional simple-cubic dis- 
persion with bandwidth W=12t, choosing i=0.5. For the 
direction of D again the y axis is selected. The impuri- 
ties have a common bath, yet V is assumed here to be k- 
independent and the explicit impurity-impurity distance 
is formally set to zero. A constant value of V~— 0.5 is 
chosen in the present work. Hence only the local part 
of the RKKY interaction is accessible. Such a model- 
ing is e.g. important for understanding the local spin 
interactions between correlated atoms on metallic sur- 
faces,— where there is indeed an intriguing interplay be- 
tween conventional direct exchange, RKKY interaction, 
Kondo effect and anisotropic exchange. 

For the numerical solution of the model Hamiltoni- 
ans discussed here, the rotationally invariant slave-boson 
(RISB) formalis m 21 ! 32 in the saddle-point approxima- 
tion, similar to the generalized Gutzwillcr approach^ 
is employed. The RISB methodology amounts to a de- 
composition of an electron operator a va with generic or- 
bital/site index fi via a jUa .=i?[^>]^/ (U v into its quasipar- 
ticle (QP) part f^ a and the remaining high-energy exci- 
tations carried by the set of slave bosons {<pAn}- Here 
A denotes a chosen localized basis state and n relates 
to the given QP degree of freedom. Two constraints, 
the first enforcing the normalization of the bosonic con- 
tent and the second keeping an eye on the match of 
the bosonic and the fermionic occupation matrix, are 
established on site-average at saddle-point through the 
Lagrange-multiplier matrix In order to describe 

inter-atomic correlations adequately, a two-site (cellular- 
cluster) framework is used. This cluster connects two 
NN lattice sites between the layers in the Hubbard bi- 
layer and the two impurities in the TIAM. It amounts 
to a local cluster approach to the electronic self-energy, 
whereby Si2(w) incorporates terms linear in frequency 
as well as static renormalizations^ Therewith the low- 
energy behavior may be adequately expressed and inter- 
site correlation functions as well as multiplet weights on 
the cluster can be retrieved. Importantly, the formal- 
ism allows for full spin and orbital rotational invariance, 
needed to account for the competition between isotropic 
and anisotropic interactions. In this respect the slave 
bosons may become true complex numbers and A can 
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be expanded via Pauli matrices in each orbital sector 
(with allowed off-diagonal terms between these sectors). 
Albeit the calculations arc formally performed at tem- 
perature T=0, a small gaussian smearing for the k-point 
integration introduces a minor T scale. For this reason 
the energetics are discussed in terms of the free energy F. 
Note that in the numerical solution of the TIAM, a three- 
orbital model is effectively treated within RISB, whereby 
the bath enters through its band dispersion. Thus the 
bath degrees of freedom are not integrated out, but are 
handled explicitly. In principle, a correlated-bath sce- 
nario may also be studied, however we here always keep 
£4>ath=0. Nevertheless, correlation effects are introduced 
within the bath due to the coupling to the correlated 
impurities. The investigated half-filled scenario of the 
model is either achieved by setting e<i=—U/2 or through 
an additional Lagrange multiplier fixing the electron oc- 
cupation on the bath according to the total filling N=3. 

III. HUBBARD BILAYER MODEL 

The original Hubbard bilayer without DM interaction 
has already been addressed in several works,— ~— most 
often concerning the electronic phase diagram when vary- 
ing the ratio t±/t. Here however the main interest lies 
on the ratio f//|D| for the coupled square-lattice layers 
with bandwidth W=8t. In the following, we restrict the 
discussion to cases t±/t<l with all the energies given in 
units of the half bandwidth At. 

Concerning the electronic phases studied within the 
current mean-field approach, wc restrict the discussion 
to local cluster orderings, i.e., neglect long-range order 
parameters suitable for e.g. spin spirals. Such more in- 
tricate instabilities are planned to be addressed in more 
concrete materials-connected future modelings. Here the 
focus is first on the interplay of the fundamental short- 
range correlation processes in the strongly correlated 
metallic regime. Note however that in the present con- 
text the cluster description does not account for intra- 
layer inter-site self-energies. 

A. Half-filled case 

At half filling, each layer accomodates one electron and 
the whole system is therefore susceptible to a Mott tran- 
sition. We study two cases, namely the one of weakly 
coupled layers (t± =0.025) and the other with stronger 
inter- layer hopping (t±=0.1). 

Figure [2] shows the phase competition within the half- 
filled model with increasing the Hubbard U. The compu- 
tations allow for the stabilization of two metallic phases, 
namely the paramagnetic (PM) and antiferromagnetic- 
between-layers (AFM) ones. From the inspection of 
the free-energy differences it is nonsurprisingly seen that 
in general the AFM phase wins over the PM phase at 
larger U. Thereby a smaller t±, and hence a smaller 



D = 0.0 D = 0.03 




FIG. 2. (Color online) Energetics, QP weight and multiplet 
weights with increasing U for the bilayer model at half filling 
for two values of t± and D, respectively. Free energies are 
normalized to the one of the PM phase. The thick lines in 
the multiplet- weights plot correspond to the states in the two- 
particle sector. 




FIG. 3. (Color online) Evolution of the Fock-state contribu- 
tions to the local state with the highest slave-boson amplitude 
(compare with Fig. [2]) for (7=3. Note that for D^O this state 
" 1" is no longer an eigenstate of the 5* 2 and the 5* z operator. 
The inset exhibits the D dependence of the state. 



bonding/anti-bonding splitting, supports the building-up 
of the AFM phase, in line with DMFT calculations em- 
ploying quantum Monte-Carlo solvers for the impurity 
problem i 20 ' 23 A further gain in AFM free energy is ob- 
served at fixed U when introducing the DM interaction, 
but with only marginal shifts of the phase onset towards 
smaller U. The difference between the two critical U =U C 
for the two different t± vanishes with D, while in the case 
of D=Q the U c for i^=0.1 is clearly smaller. However the 
general evolution of the QP weight Z=(l— 9£/i9u;) -1 | u= o 
with U does not display strong changes with the intro- 
duced anisotropic interaction. 
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In addition, Fig.[3]also displays the slave-boson weight 
of the local multiplet state that dominates at half filling 
in the relevant two-particle sector. For D=0 the bilayer 
Hamiltonian (JlJ commutes with {S 2 , S z } and thus singlet 
and triplet states form the local two-particle eigenbasis. 
As expected, the singlet remains strongest up to the Mott 
transition, followed by the triplet states (whith their de- 
generacy lifted when entering the AFM phase). In the 
case of finite D, the picture formally looks very similar, 
but Hamiltonian and spin operators are no longer com- 
muting operators and the respective two-particle states 
thus are not true spin eigenfunctions. That is easily un- 
derstood from the DM interaction favoring a perpendic- 
ular alignment of the local spins, contrary to the origi- 
nally preferred collinear states. It is nicely illustrated in 
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FIG. 4. (Color online) Interaction dependence of the spin 
expectation values in the AFM phase. Nonzero values of the 
DM coupling introduce a spin component pointing along the 
x axis (due to the choice for the direction of the D vector). 
The vertical dot-dashed lines mark the AFM transition point. 
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FIG. 5. (Color online) Same as Fig. [4] here for the spin- 
spin correlation functions between the layers. Solid lines: PM 
phase, dashed lines: AFM phase. 
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FIG. 6. (Color online) Influence of the DM coupling on the 
QP weights, spin moments and spin-spin correlation functions 
for the AFM half-filled bilayer model. The inset shows the 
evolution of the angle between the spin moments in the two 
layers. 



Fig. [3J where the inter-mixing of the singlet and triplet 
Fock-state building blocks with growing D is exhibited. 

Figure U shows the evolution of the spin moments in 
the two layers with increasing U. For D—0 only (S z ) 
adopts a nonzero value in the AFM phase, with a steeper 
increase for larger t±. However with finite D also a size- 
able x component of (S) shows up and grows until U c 
is reached. For the smaller t^=0.025 the value for (S x ) 
even equals the corresponding (S z ) magnitude. A lower 
t± apparently also effectively increases the relative ten- 
dency towards the non-collinear spin alignment driven by 
the DM coupling. Note that the DM interaction not only 
modifies the AFM phase, but has an impact in the PM 
state as well. There (Si^S^y) exhibits more AFM-like 
character and the corresponding (x, z) correlation func- 
tions show minor weakened AFM-like tendencies, both 
compared to the D=0 case. Close to the Mott transi- 
tion the larger t± results in a stronger (coherent) spin 
response for D=0, as retrieved from the inter-layer spin 
correlation functions plotted in Fig. [5] For nonzero D the 
correlation between the x components, i.e. (Si^S^x), 
appears to behave especially more disconnected from the 
z component for the smaller t±. 

In order to gain further insight into the impact of the 
DM term, Fig. [5] depicts explicitly the D dependence for 
fixed U. The Mott transition itself may be tuned over a 
rather wide range of the anisotropic interaction. Whereas 
the spin moment in the x direction shows a strong vari- 
ation with D, the spin-spin correlations are only weakly 
dependent thereon. Albeit no resulting (S y ) value exists, 
the correlations along y still appear to gain strongest in 
magnitude from a larger D. It is also visualized that 
the angle 7 between the local spins on the adjacent lay- 
ers indeed increasingly deviates from the AFM-ideal 180° 
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with growing DM interaction. Close to the Mott transi- 
tion, the value 7^100° is nearby the DM-ideal value of 
90°. 



B. Hole-doped case 

We now turn to the effects of doping the bilayer model 
away from half filling. For investigating the simultaneous 
effects of doping, on-site Coulomb interaction and inter- 
site DM interaction we set t±—0.1 and first fix the Hub- 
bard interaction to f/=3. As can be seen from Figf2]thc 
value of U puts the system just below the Mott transition 
at half filling, i.e. strong correlations with the quasipar- 
ticle weight exist. 

The results of hole doping 8=2— n for the system in the 
filling range n G [1.6, 2.0] are summarized in Fig. [7] Let 
us first discuss the D=0 case. Starting from half filling, 
the system is in the AFM phase for the chosen U value. 
With increasing 5 the order parameter (S z ) decreases, un- 
til it vanishes close to n=1.74 and the PM phase sets in 
(at reduced spin-spin correlations and larger QP weight). 
When including a DM interaction with Z?=0.03 in the 
model, the situation becomes more intriguing. Again the 
AFM (C-AFM) phase, now canted in x direction, weak- 
ens upon doping from half-filling, however at n^l.76 the 
Hubbard bilayer system shows a first-order phase tran- 
sition to a metallic spin-flop (SF) phase. The latter 
one is characterized by the discontinuous jump to a lo- 
cal configuration with an (S x ) expectation value larger 
than (S z ). This corresponds to an angle 7 between the 
local spins in both layers being lower than 90° , whereas 



D = 0.0 D = 0.03 
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FIG. 7. (Color online) Influence of hole doping on the bi- 
layer model, with and without DM interaction for t±=0.1 
and U=3. Insets at the top show the free-energy curve, with 
the region where AFM order is (meta)stable marked in red. 
C-AFM marks the canted antiferromagnetic and SF the spin- 
flop phase. The inset in the right-middle depicts the evolution 
of the angle 7 between spins. 




PM AFM C-AFM SF 

FIG. 8. (Color online) Illustration of the stable local spin 
configurations according to the local spin-spin angle 7 on the 
inter-layer cluster, (a) PM without ordered local moments, 
(b) AFM with 7=180°, (c) C-AFM with 7 between 180° (pure 
AFM ordering) and 90° (pure DM ordering), (d) SF with 
7<90°, i.e., weak ferromagnetism with strong canting. 

in the C-AFM phase 7 G [90°, 180°] holds (see Fig©. 
The strong decrease of 7 at the transition point may be 
observed in the inset of Fig. [7J Hence the SF phase dis- 
plays weak ferromagnetism due to strong canting. Note 
that neither the spin correlation functions nor the diag- 
onal Z values show a strong signature therein. The SF 
phase transforms into the usual PM phase at n~1.62. 

In addition to the doping scan, Fig. [5] displays the var- 
ious phases emerging with increasing Hubbard interac- 
tion U for fixed hole doping (5=0.3, i.e., n=1.7. Without 
the DM interaction, the standard picture of a stable PM 
phase at small U and a stable AFM phase at larger U 
(£/>3.58) remains vital. Note that the U values for AFM 
stabilization are well above the Mott critical U at half fill- 
ing. Introducing D stabilizes the metallic SF phase for 
1.9<C/<3.85, accompanied with the jump in the angle 7 
towards lower values. Therewith the onset of AFM order 




u 



FIG. 9. (Color online) Phase diagram with U for the doped 
Hubbard bilayer at filling n=1.7 (compare with Fig. [7]). In- 
sets show again the free-energy curves and the interaction- 
dependent spin-spin angle 7. 
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takes place at slightly larger U than for D=0. Hence the 
finite D enables specific magnetic ordering in a Coulomb 
interacting regime that is originally not susceptible to 
such order. Only the z component of the spin correlation 
function shows a discontinuous behavior at the SF/C- 
AFM phase boundary. 



IV. TWO-IMPURITY ANDERSON MODEL 

The TIAM2£~— belongs to the set of canonical models 
in the physics of strong electronic correlations, believed 
to be relevant for the understanding of heavy- fermion sys- 
tems^. Via the coupling of the impurities to a bath it 
contains the single-impurity Kondo physics and as com- 
petitor also the RKKY mechanism acting between the 
impurities. The latter originates from the effective ex- 
change introduced through the impurity-coupling to the 
same bath. In some work a 18 ' 31 this type of exchange 
interaction between sites is discussed in the context of 
two impurities coupled to different baths (similar to the 
bilayer architecture). But here we try to separate the 
exchange in an indirect ("RKKY") one, stemming from 
effective exchange via the bath, and a direct term, result- 
ing e.g. from an explicit hopping amplitude t\i between 
the impurities (see Eq. ([2])). 

Figure [TU] shows the RISB results for the spin correla- 
tion functions of the fundamental model with t\2=D=0. 
Thus the two impurities arc only coupled via the bath 
and exchange can only be mediated therewith. The ex- 
pectation value (S 2 )=S(S + 1) on the impurity quickly 
rises with U due to the formation of the local moment. 
It approaches the value 3/4, corresponding to the full 
5=1/2 limit, at large interaction strength. With increas- 




FIG. 10. (Color online) Interaction-dependent spin correla- 
tion functions for the TIAM with ti 2 =0. Top: (S 2 ) for D=0, 
bottom: spin-spin correlations for D=0 (circles) and D=0.05 
(squares). The dashed dark line without circles is just the 
(S'impS'bath) correlation function mirrored at the zero line. 




FIG. 11. (Color online) Interaction-dependent spin-spin cor- 
relations for the TIAM with t 12 =0.2, with and without DM 
interaction. 



ing U a local Fermi liquid is established with a small 
quasiparticle weight Z- lmp (see inset Fig. [TU|) . The compe- 
tition between the Kondo screening and the RKKY inter- 
action may be observed from inspection of the spin-spin 
correlations. From Fig.[TU]it is obvious that (S'impS'bath)! 
i.e. the correlation between a single impurity and the 
bath, is always of AFM character with a maximum close 
to Z7k~1-6. On the other hand the inter-impurity corre- 
lation (S1S2) is exclusively of FM kind and shows mono- 
tonic increase with U . The former is associated with 
the singlet-forming tendencies due to Kondo screening, 
whereas the latter signals triplet-forming tendencies be- 
cause of the FM RKKY exchange within the local limit. 
Close to C/k the absolute value of the local RKKY cor- 
relation exceeds the singlet-forming amplitude between 
impurity and bath. The system at larger U is then dom- 
inated by the RKKY interaction.^ .^ Within a conven- 
tional Schrieffer- Wolff mapping 3 ^ for the Kondo coupling 
via Jk=8V 2 /U , a similar crossover regime would follow 
also from simple estimates through the associated ex- 
change interactions. For if we understand the RKKY in- 
teraction as second-order process, i.e. Jrkky~^kj then 
here the two exchange integrals become equivalent for 
[7=2, which is the order of magnitude from the nu- 
merics. With increasing impurity-bath coupling V the 
crossover shifts to larger U, since Jj| stronger profits 
therefrom. However note that with our bath bandwidth 
IF^G the present TIAM is surely not in the Kondo- 
Hamiltonian limit (U^>W) for the studied interaction 
ranged Turning on a finite DM term of size -D=0.05 has 
nearly no effect at small U . However for larger Hubbard 
interaction rather strong modifications occur especially 
for the inter-dimer function (Si, y S2, y ). Remember that 
the D vector also points in the y direction. Thus an in- 
tricate spin-spin coupling scenario arises at large U , with 
still FM alignment in the (x, z) axes and near AFM align- 
ment in the y axes. For U>5 our mean-field approach 
yields net local moments in presence of a finite D, i.e., a 
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paramagnetic solution is no longer stabilizable. It would 
thus be very interesting to study the large- U regime of 
this model beyond mean field (e.g. with the numerical 
renormalization group approach utilized in Ref . l30l) . 

In addition to the basic model with vanishing inter- 
impurity hopping, Fig. [TT1 exhibits the resulting spin cor- 
relation functions for the TIAM with ti2=0.2. Now both 
(S'impS'bath) an d (S1S2) display AFM correlations in the 
weakly interacting limit. This is understood from the di- 
rect exchange integral Jdh=4t\ 2 /U originating from the 
introduced dimer coupling. With increasing U the corre- 
lation functions develop rather similar as for t±2=0, yet 
the overall magnitude is somewhat reduced a small inter- 
action strength. Hence there the direct exchange weak- 
ens both, impurity Kondo-screening (due to the stronger 
inter- impurity link) as well as FM RKKY interaction 
(since the direct exchange favors AFM behavior). But 
the crossover point of domination for these processes does 
not seem to change much with the introduced ti2- Of 
course, a very large ti2 should rank the direct exchange 
above the other mechanisms, however here we do not 
investigate this model limit. Finally, when introducing 
the DM term to the model, effectively four different ex- 
change mechanisms compete with each other: impurity 
Kondo, RKKY, direct and DM. The latter has indeed 
again significant effect on the spin correlation between 
the impurities. For already moderate values of U the 
dominance of the FM RKKY is lost, turning the sys- 
tem into AFM-likc inter-impurity correlations for U>2. 
Thus also here the DM interaction severely influences 
the magnetic correlations for isolated impurities within 
an itinerant background. It appears to strengthen the 
singlet-forming tendencies (with stronger response in the 
D direction) in an otherwise triplet-favoring RKKY sys- 
tem at short-range distance. 

V. SUMMARY 

A theoretical investigation of effects stemming from 
the Dzyaloshinskii-Moriya interaction in itinerant sys- 
tems with strong electronic correlations was presented 
in this work. In order to study the principle physics on 
the lattice as well as in the local limit, we elaborated 
on two prominent model systems, namely the Hubbard 
bilayer and the one defined by the two-impurity Ander- 
son Hamiltonian. In both cases substantial influence of 
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the DM interaction is found, especially at large coupling 
where the strong renormalization enhances the impact. 
The half-filled Hubbard bilayer displays large out-of-axis 
spin components close to the Mott transition that may 
severely influence the magnetic response in applied field. 
Intriguing phenomena in this respect are e.g. observed 
in quasi-two-dimensional organic compounds^ At finite 
hole doping and larger U, the bilayer system with DM 
interaction exhibits the emergence of a metallic spin-flop 
phase inbetween the AFM phase at half filling and the 
PM phase at stronger doping. This finding is of vital im- 
portance for many doped Mott systems with anisotropics. 
For instance, it is well-known that the hole-doped layered 
cupratcs display puzzling phases inbetween the AFM and 
the superconducting dome and that the DM interaction 
is not completely negligible at low energy^ - — Thus it 
would be very interesting to investigate in some detail 
whether there is a closer connection between our model 
results and those observations. 

The TIAM is very relevant not only in the context of 
heavy fermions, but e.g. also for isolated atoms on metal- 
lic surfaces. In either case, anisotropic spin terms such 
as the DM interaction exist in many realistic represen- 
tants in nature. It results from our studies that the DM 
term becomes an important player in the hierachy of rel- 
evant exchange processes in these contexts. In the local 
limit it works against the FM tendencies of the RKKY 
interaction and promotes the singlet formation between 
the impurities at large local Coloumb interactions. Fur- 
ther research along these lines, e.g. by going beyond 
mean-field, including the complete k dependence of the 
impurity-bath coupling or tailoring the modeling towards 
concrete materials systems is of vital interest to account 
for generic exchange processes in the strongly correlated 
itinerant regime. 
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